Solvable model for the standard folklore of the glassy state 
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A model system with last and slow processes is introduced. Alter integrating out the last ones, 
the considered dynamics ol the slow variables is exactly solvable. In statics the system undergoes 
a Kauzmann transition to a glassy state. The relaxation time obeys a generalized Vogel-Fulcher- 
Tammann-Hesse law. The aging dynamics on the approach to and below the Kauzmann temperature 
is analyzed; it has logarithmic behavior. The structure ol the results could be general, as they satisly 
laws ol thermodynamics lar from equilibrium. The original VFTH law is on the border-line between 
the regime where only the effective temperature ol the slow modes is needed, and the regime where 
also an effective field occurs. The production rates ol entropy and heat are calculated. 
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Half a century ago Kauzmann pointed at the "entropy 
crisis" in glass forming liquids: The naive extrapolation 
of the entropy of an undercooled liquid goes below the 
entropy of the crystalline state, which is physically incor- 
rect M . It is often assumed that a thermodynamic phase 
transition takes place around the "Kauzmann tempera- 
ture" Tk-i characterized by a vanishing configurational 
entropy. This prediction is very difficult to test exper- 
imentally, since the relaxation time may be very long. 
Experimental data for glass forming liquids are often fit 
to a Vogel-Fulcher-Tammann-Hesse (VFTH) || behav- 
ior T cq ~ exp[v4 7 /(T — To) 7 ]. An entropic argument for 
the standard choice, 7 = 1, was given by Adam and 
Gibbs H . A field-theoretic approach by Parisi ^ leads to 
7 > 1. These approaches share the commonly accepted 
point that To and Tk coincide. In practice exponents 
7 7^ 1 are compatible with data; this merely affects the 
width of the fitting interval. Here we shall investigate 
universal aspects of this standard picture. 

Very little is known about dynamics below the Kauz- 
mann temperature. One simple, intriguing idea comes 
to the mind. Consider a glassy system that has aged 
a long time t. Can one say that all processes with 
equilibration time much less than t are in equilibrium, 
while the ones with timescale much larger than t are still 
quenched, leaving the processes with timescale of order t 
as the only interesting ones? Such an approach would 
be too naive for continuous phase transitions, where 
the timescale diverges algebraically at the critical point. 
However, in glasses there is an exponential divergence, 
and this might induce an asymptotic decoupling of the 
time-decades. Indeed, it occurs in the glassy regime of 
models with an Arrhcnius (rather than VFTH) law, like 
the ferromagnetic Ising chain |5j, harmonic oscillator 
and spherical spin models ||. This asymptotic decou- 
pling could be the basic ground for a generalization of 
equilibrium thermodynamics to systems far from equi- 
librium || . That approach involved systems where one 
extra variable was needed to describe the non-equilibrium 
physics, namely the effective temperature. In the present 



work we shall encounter situations, where also an effec- 
tive field is needed to describe the two observable quanti- 
ties, energy and magnetization out of equilibrium. That 
case is of principle interest for glass forming liquids, as 
it was found experimentally that at given T, p the same 
volume can be reached via different histories. As they 
then lead to different futures, the glassy state of these 
systems cannot be coded using one extra variable Q . 

In order to gain insight in properties of the standard 
picture, we introduce a simple model, with solvable dy- 
namics. In the long time aging regime it leads to first 
order differential equations. We thus can study aging on 
the approach to and below the glass transition. Even 
though the physics of the model is quite simple, we shall 
set out to discover general aspects of the results by phras- 
ing them in a thermodynamic language. 

The model contains a set of N spins Si, coupled to 
a slowly changing background field, that we model by a 
set of N harmonic oscillators Xi. All variables are only 
coupled locally, as is described by the Hamiltonian 



H 



(1) 



The spins have no fixed length, but satisfy the spherical 
constraint JV &i = N. In order to model the /3-processes 
of a glass, we assume that the spins Si move fast, and 
we shall consider only the time regime where they are at 
thermodynamic equilibrium. To model a-processes, we 
shall consider a dynamics where the Xi may fall out of 
equilibrium. Integrating out the Si replaces the spin-part 
of the energy (Q) by its free energy, which will act as an 
effective Hamiltonian for the slowly moving variables Xi, 



T~tcS — U — TS, 



cp 



(2) 



Here U = Nu is the internal energy given by 

U l r, rr T 

u = — = —A mo — tlm\ — w H , 

N 2 2 

It depends on the two lowest moments of the Xi, 
(l/N)M k = (1/N) Y,i x t and w is the short hand 
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2JL mi + L 2 + -T 2 , 
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5 cp is the entropy of the equilibrium processes, i.e., of 
the spins, 



S cv = l -N- l -N\n[ 



(•5) 



Being interested in the long-time aging behavior, we may 
restrict ourselves to the time-regime where the dynamics 
of the Xi can be described by a Monte Carlo dynamics. 
We must assume parallel updating, as sequential updat- 
ing would not produce a long timescale in our system 
without two-particle interactions. We are thus consider- 
ing a simple example of models typically studied numer- 
ically. Since it can even be solved analytically, we shall 
gain qualitative insight in the long time regime. 

The rules for the MC dynamics are standard, and lead 
to evolution equations for mi 2 that have the form 



00 

'"1 : / dxyi(x)p(x\mi^)W{[3x) 
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m 2 =2 I dx X + Hyi ^ P (x\ mi . 2 )W(f3x) (7) 

J-oc K 

Let us explain how this comes about. For one MC 
step one shifts all Xi by an amount ri/y/W with the 
Ti drawn from a symmetric Gaussian with variance A 2 . 
This updates M{ 2 = M1.2 + SMi, 2 . According to the 
Metropolis algorithm, the move is accepted with proba- 
bility W(f3x) = exp(—(3x) if the change in effective en- 
ergy, x = H' cS — Ties-, is positive; if x is negative the 
move is always accepted (W((3x) — 1). Such dynamical 
problems are well understood || || . The essential ingre- 
dient is the thermally averaged transition probability for 
a move characterized by x and 8M\ . It is the product of 
two factors: the first is a Gaussian in 8M\ with center 
yi{x) and width A y ; the second factor, p{x\mi y2 ) occur- 
ring in eqs. (|^) and (Q), is a Gaussian in x with center 
xq and width A^. Defining \i\ = H / K — m\ with 



H = H 



JL 



K = K 



J 2 



(8) 



w + T/2 1 " " w + T/2 

the parameters of these Gaussians take the form 

x Q = A 2 K/2, A x = A 2 KT e , (9) 

2/i = /3 e /ii (x - x), A y = A 2 (1 - e Kiil). (10) 

Here T e = l/(3 e is a short hand for T e = K(m 2 — mf); 
later we shall interpret it as an effective temperature. 

We are still free to choose the variance A 2 of the Monte 
Carlo moves. To model glassy behavior we shall consider 
a dynamics that constrains the system's phase space to 
m 2 — m\ > mo for some fixed mo- This mimics a glass 
that has no crystalline state, as occurs in atactic polymer 
melts. For some fixed mo and B, we thus let A depend 
on the instantaneous values of mi^t) according to 



A 2 (t) = S^r(t) 

V ; K{t) V ; 



(11) 



where 

r(t) = 



( B 



Ma(*) = m 2 {t) - m\{t) - m (12) 



The power 7 remains a free parameter, because we will 
not make a connection to an underlying microscopy. 

The fact that A 2 becomes large near \i 2 = implies 
that there the trial moves are typically large, and thus 
unfavorable. Most of them will be rejected, implying, as 
shown below, a VFTH-type equilibrium relaxation time 



T eq (T) ~ expT e 



exp 



( A(T ) 

It- T 



(13) 



while T e „ — 00 for T < Tq. Though the physics of our 
model (0) is simple, the choice (pT|), ( |l2| ) for the timescale 
allows us to derive aspects of the long time aging regime 
below Tk, which could well be model independent. 

Let us first look at the static regime, that has \i\ = 0. 
For large enough T, the variable fi 2 remains positive and 
the effective temperature is trivial: T£ q = T. A Kauz- 
mann transition occurs at Tr- set by /x 2 = 0. Since states 
with \x 2 < cannot be reached, it is evident that the 
generally expected equality To = Tk holds. In statics [i 2 
vanishes for all T < Tq, implying that T° q follows from 



Km - T e oc i + T D 2 m 



J 2 m 



KniQ — T, 
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( JT e eq ) 2 (Km - T e eq ) 2 
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where D = HJ + KL. The Kauzmann temperature 
is the solution of T^(T ) = T . Below T , T^ ex- 
ceeds T, and it remains finite for T — > 0. All static 
observables can now be expressed in T° q . Since w = 
[J 2 m + (DJ/TfJ) 2 +T 2 /4] 1 / 2 , we have the result for 
mi = H/K via (@), and using m 2 — m\ + mo, we can 
find the energy from (§). Since dT^/dT < 1 at T ", the 
specific heat makes an upward jump as function of T; 
this happens also for —dmi/dT. The very same features 
occur in the mean field glass model, the p-spin model, 
where T^ = x{T)/T, with x the Parisi breakpoint (k]]. 

We shall restrict ourselves to aging experiments at 
fixed T, H, where the system evolves from some initial 
state. An example is fast quenching from an equilibrium 
state at large temperature, for which initially [i\ = 0, 
< fx 2 <C 1. In the long time regime eqs. (||), (0) lead 
to first order differential equations for fj,i l2 



fix = 4[aJrT e - (1 + a£>)(r + r)/xi]7 
rT 



(15) 
(16) 



where we denote r=(T e - T)/(2T e 
and we introduce the shorthands 



T), I = / dxpW, 



2 



DJ 2 



K 3 w(w + T/2f 



b = 



J 4 T e 



2K 2 w(w + T/2f 



(17) 



fi2 in eq. ( |16[ ) is dominated by the small factor / ~ 
exp(— r). This is the dynamic imprint of the static VFTH 
law (|l^), which may be more general. The solution fol- 
lows readily, since the integral of an exponential is basi- 
cally the same exponential. It is this property that leads 
to the asymptotic decoupling of decades, discussed in the 
introduction. In the aging regime one gets t ~ to exp T, 
equivalent to 



D 



M2 
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where to depends logarithmically on t. Such logarithmic 
behavior is supposed to be the universal finger print for 
the glassy state; notice, however, that there can be a non- 
trivial exponent. For T > Tq fi2 will stick at ^ ~ T—Tq, 
whereas for T < Tq it will ultimately vanish. 

To solve for /Ui(t) we divide eq. dl5| ) by (|l6|), which 
allows to express \i\ in terms of // 2 of eq. (pi 



d/ii 
d^2 



-aJKrT e + (1 + aD)(T + r)Kfj,i 
2rT e - 2(r + r)it(4 



(19) 



There exist several ranges of solutions: 
• Regime 0: T > T , all 7: 

Here r ~ T e — T ~ ^2- The solution is obtained by 
equating the numerator of the right hand side to zero, 



Mi 



aJK(l 



1 + a£> T 



(20) 



Since /xi <C ti 2 , it can be neglected to leading order. 
• Regime 1: T < T , 7 > 1: 
As in Regime 0, but now r is of order unity, 



/'1 



a JT e r 
1 + aD F 



(21) 



Also in this case n\ <C /i2- 

• Regime 2: T < T , 7 = 1: 

This is the case with a true VTFH-law. Quite surpris- 
ingly, eq. ([l9]) exhibits several subregimes with different 
behavior: 7 = 1 is the most difficult situation, too intri- 
cate to explain here in detail. 

• Regime 3: T < T , 7 < 1 

Now the denominator of ( |l9| ) becomes subleading, of 
order T/i2 = B 1 pb 2 7 <C 1, but it remains positive, thus 
assuring that \ii finally vanishes. It follows that 



so in this regime /x 2 <C /ii can be neglected. 
The magnetization approaches equilibrium as 



H cq 1 



aJK 



2(1 + aD) 



Ma(*) (23) 



In regimes and 1 its time-dependence is dominated by 
/X2 ~ (hit) -1 / 7 , while in regime 3 the dominant term is 
[i\ ~ (lnt) -1 / 2 . In all cases both m,2 = mo + m\ + /i 2 
and the energy, see eq. (^), behave similar to m\. 

It has long been a challenge to formulate glassy behav- 
ior in terms of the laws of thermodynamics, because that 
could point at universal behavior. Attempts to do so re- 
mained unfruitful till recently |j |ll]] . In a recent series 
of papers we have shown that it can be done in model 
systems where only one extra variable, the effective tem- 
perature T e , is needed to describe the physics [pjj || 
||. The crucial step was to explain the paradoxes con- 
cerning the Ehrenfest relations and the Prigogine-Defay 
ratio pl| . In the present model there are two indepen- 
dent observables, mi^ft) or, equivalently, m\{t) and u{t). 
They result from the dynamics, but they can also be de- 
scribed by a thermodynamics at an effective temperature 
T e (t) and an effective field H e (t). To do so, let us cal- 
culate in a quasi-stationary approach the partition sum 
over macroscopically equivalent states 

Z e = [ Dxe-^ H °«+ HM i- H ° M ^S(M lt 2 - Nm h2 (t)) 



(24) 



By the extra terms in the exponent we eliminate H in 
favor of H e . We find for the free energy F = —T e \nZ e 



F = U + (H- H e )Mi - TS cp - T e l 



(25) 



where U and S'op were defined in (§) and (§. The new 
term I = N[l + ln(rn2 — mf)]/2 is the entropy of the 
configurations with given m\ t 2- Though our model is 
classical, we can impose X = at T — by adding a 
constant, yielding I — (N/2) ln(l + ^/ttiq). In his sem- 
inal paper Kauzmann pointed out that the equilibrium 
configurational entropy vanishes already at Tk- We see 
that this is indeed the case here. We can also check an- 
other aspect of the standard folklore. Comparison with 
eq. ( p"3| ) yields lnr eg ~ (N/J) 1 . This is the generalization 
of the Adam-Gibbs relation jj| to cases with 7 7^ 1. 

T e and H e are required to satisfy the saddle point equa- 
tions that follow from varying m\ : 2 in eq. (pa), implying 



T e = k(m 2 - ml) = K(m + (j, 2 ) 



(26) 



H e = Km 1 - J -^± J ^ = H-k, l (27) 

They only depend on H through m\ t 2 — mi r2 (t;T, H). 

In eqs. (p|)-(|TT|) it was already seen that the expression 
( ^6| ) for T e was useful to shorten the dynamical equa- 
tions. It can now be checked that T e — T° q ~ +fi 2 — Mi > 
so it has no definite sign, whereas H e — H ~ — /ii is neg- 
ative. It is now seen that the "simple" Regimes and 
1 have a natural interpretation: one then has \i\ <C /i2j 
so the effective field H e w H is not needed: the glassy 
state is coded by the effective temperature alone. Such 



3 



a situation was encountered before in model glasses with 
T = § §, and in the p-spin model @ @. In 
the "complicated" Regimes 2 and 3 one ultimately has 
T e — T£ q ~ H e — H . Both parameters are thus needed to 
describe the evolution, a result that should be valid be- 
yond our model. Systems with a true VFTH-law (7 = 1) 
are at the borderline, and thus more subtle. 

We now consider which processes are involved at large 
time t. The change of energy may be cast in the form 



u(t) 



u(t) 



-(T,T e (t)) 



(28) 



For T > Tq the behavior I ~ exp(— T) reproduces at 
equilibrium eq. (|l3|); Outside equilibrium it follows that 
t(T, T e ) w r eq {Te), as was observed in systems with an 
Arrhenius law gj. For T < T , 7 > 1 (Regime 1) the 
timescale of a system characterized by T and T e obeys a 
generalization of the VFTH-law 



t(T, T e ) ~ exp r = exp 



A(T) 



r e - T e cq (T) 



(29) 



saying that the equilibrium VFTH law has a definite fin- 
gerprint on the aging dynamics below To. Here, as well 
asineq. ©, A{T) = BK(l + 6/(1 + aD)). In the cases 
7 < 1 there are various behaviors, that may not be uni- 
versal. In all these cases it holds that for t 3> to the active 
processes are indeed those with timescale r(T, T e (t)) ~ t, 
as is consistent with the logarithmic time dependence: 
u(t) = u(lnt) implies u = —u/t — ► r = {—u/u')t. 

Let us now consider the results from a more general 
thermodynamic point of view. The first law says that the 
change in energy is equal to the heat added to the system 
and the work done on the system, dU = &Q — Midi?. It 
is satisfied with 



dQ = TdScp + T e dl + (H e - H)dM x 



(30) 



This holds because of the presence of a free energy, eq. 
@, satisfying dF = -S cp dT - ZdT e - MtdH e . It does 
not invoke the relation between m\{t) and m2(<), which 
differs in the various regimes. The two entropic terms 
in ( ^o|) are related to the fast and slow processes, re- 
spectively B p^] . The last term is new; it results from 
irreversible internal work. 

In the aging process the system will go to lower energy 
by producing heat at a rate |dQ/dt|, that follows from 
©. It holds that |dQ/dt| ~ |/ii| + \fi 2 \ • The dominant 



term is S, 
due to S, 



ep 



op 



- X ~ fi2 in regimes and 1 , while it is solely 
' /ii in regime 3. 
In thermodynamics close to equilibrium a central role 
is played by the irreversible entropy production. For sys- 
tems far from equilibrium, this is also the case. We now 
show that it can be derived explicitly for systems such 
as the present one. We write the change in total entropy 
S = S cp +2 as dS = d e S + diS, where d e S = dQ/T is 
the externally supplied entropy, and diS the internally 
produced entropy. This gives the general result 



dSi 
dt 



T-T e dl 
~dl 



H -H e dMi 



(31) 



T dt T dt 

It can be verified explicitly in our model that both contri- 
butions are non-negative, in accordance with the second 
law. The first term of (|3l]) is of order rfi-z, and is domi- 
nant in regimes and 1; The last term is of order ^tx|Ai|) 
and it dominates in regime 3. 

In conclusion, we introduce an exactly solvable model 
glass with a set of fast variables (spherical spins) and 
a set of slow variables (harmonic oscillators). The fast 
modes are summed out. For the slow modes we as- 
sume a dynamics that obeys a generalized VFTH law, 
r eq ~ exp A~*(T — Xo)~ 7 . In statics there occurs a Kauz- 
mann transition. We solve the aging dynamics, and ver- 
ify a generalized Adam-Gibbs relation between the re- 
laxation time and the configurational entropy. Below To 
the relaxation time for a system described by its effective 
temperature may have a VFTH form, or a different one. 
The dynamically active processes have timescale ~ t. 

The aging solution can be described by a quasi-static 
thermodynamic approach. Above the Kauzmann tem- 
perature To one extra variable is needed, the effective 
temperature, that depends logarithmically on time. Be- 
low To this persists in some cases, while in other cases 
also an effective field is needed. Systems with the original 
VFTH-law are at the borderline. 

We finally derive the production rates of entropy and 
heat for systems described by an effective temperature 
and field. They are non- vanishing whenever the effective 
temperature deviates from the bath temperature, and/or 
when the effective field deviates from the external field. 
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